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ABSTRACT 

We elaborate on a recently proposed model for subsonic quasi-spherical accretion onto slowly 
rotating pulsars, in which accretion is mediated through a hot quasi-static shell above the 
neutron star magnetosphere. We show that under the same external conditions, two regimes 
of subsonic accretion are possible, depending on if plasma cooling in the transition zone is 
dominated by Compton or radiative processes. We suggest that a transition from the higher 
luminosity Compton cooling regime to the lower luminosity radiative cooling regime can be 
responsible for the onset of the 'off' -states repeatedly observed in several low luminosity 
slowly accreting pulsars, such as Vela X-l, GX 301-2 and 4U 1907+09. We further suggest 
that the triggering of the transition may be due to a switch in the X-ray beam pattern in 
response to a change in the optical depth in the accretion column with changing luminosity. 
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1 INTRODUCTION 

Quasi-spherical accretion from stellar winds onto magnetized neu- 
tron stars in binary systems may proceed in two physically distinct 
ways. The gravitationally captured stellar wind matter is heated be- 
hind the bow shock at the characteristic Bondi radius Rg ~ 2GMjv\ v 
where v w is the relative wind velocity. If this matter cools down on a 
time scale shorter than the free-fall time tff = R 3/2 / V2GM, it will 
fall at a supersonic velocity towards the magnetosphere and come 
to a halt in a shock above it. The plasma cools down most effec- 
tively via Compton scattering of X-ray photons generated near the 
neutron star surface. Therefore for such supersonic (Bondi-type) 
accretion to take place a dense photon field is required. Thus, this 
regime of accretion is supposed to take place at high X-ray lumi- 
nositie s and was stud i ed in detail, for example, in lArons & Leal 
(1976): Burnard ct al. I dl983h . In the Bondi accretion regime the 
matter cannot pile up above the magnetosphere, and the mass ac- 
cretion rate onto the neutron star M is ultimately determined by the 
Bondi-Hoyle-Littleton formula for the gravitational capture mass 
rate by a moving neutron star M m p w R 2 B v K . 

At low or moderate X-ray luminosities the captured wind mat- 
ter may have no time to cool down, the mass fall rate becomes sub- 
sonic, and a hot quasi-spheri cal shell is due to be form ed above the 
neutron star magnetosphere (Davies & Prinddfl98lh . At the base 
of t he shell the plasma m ust cool down to some critical temper- 
ate ( Eisner & Lambll 19771) in order to enter the magnetosphere via 
instabilitites, and the rate of plasma entry into the magnetosphere 
will be determined by these instabilities. Therefore, the velocity 
of matter settling through the shell will be regulated by the abil- 
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ity of the plasma to enter the magnetosphere. As was shown in 
( Shak ura et alj|2012l Paper I hereafter), the extended quasi-static 
shell mediates the angular momentum removal from the rotating 
magnetosphere by large-scale convective motions. The mass accre- 
tion rate onto the neutron star is determined by the density above 
the magnetosphere and the mean settling velocity of matter through 
the shell, and can be very small if plasma cooling above the mag- 
netosphere is inefficient. 

It was shown in Paper I that a settling regime of accretion onto 
slowly rotating neutron stars can be established for X-ray luminosi- 
ties L x < L* as 4 x 10 36 erg/s corresponding to accretion rates 
M < M* = 4 x 10 16 g/s. At higher accretion rates, a free-fall gap 
above the magnetosphere appears in the flow due to rapid Compton 
cooling, and accretion becomes highly non-stationary. 

The model of subsonic quasi-spherical settling accretion pre- 
sented in Paper I is based on generic properties of wind accretion 
onto magnetized neutron stars with moderate and low X-ray lu- 
minosities, and can be applied to a variety of sources. For exam- 
ple, in Paper I we applied this model to observations of the slowly 
rotating X-ray pulsars Vela X-l and GX 301-2 in, both in high- 
mass X-ray binaries and spinning at an equilibrium period, as well 
as to the steadily spinning-down X-ray pulsar GX 1+4 in a sym- 
biotic X-ray binary, in which ne gative spin-down-luminosity cor- 
relations are observed (see also iGonzalez- Galan et al 112012b . The 



model was also successfully used by Lutovinov et alj 1 20 1 2) to de- 
scribe the observed spin-luminosity correlations in the slowly ro- 
tating low-luminosity X-ray pulsar X Per in an Be/X-ray binary 
as we ll as the observed lon g period in the Be/X-ray binary SXP 
1062 jPopov & Turollal2012!) . The model was also recently used in 
popul ation synthesis studies dChashkina & Popovl 120121 ; Lit ct al. 
120121) . 
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In this paper we further develop the model by focusing on two 
different regimes (within our subsonic model) of plasma entering 
the neutron star magnetosphere. In Section 2 we show that under 
the same external conditions in the wind (p w , v„ ), which largely de- 
termine the plasma density distribution in the quasi-spherical hot 
shell, subsonic accretion through the shell can occur in two dis- 
tinct regimes depending on the characteristic cooling time of the 
plasma above the Alfven surface: the Compton (shorter time scale, 
higher luminosity) and radiative (longer time scale, lower luminos- 
ity) regime. The plasma cooling time t coo \ determines the mean ve- 
locity of matter falling through the transition zone. In this zone, 
the Rayleigh-Taylor instability, which allows plasma to enter the 
magnetosphere, develops. This velocity is inversely proportional to 
the plasma cooling time, u R ~ f Q , 3 , and eventually determines the 
mass accretion rate through the magnetosphere M onto the neutron 
star. Therefore, in the same source accretion in both a higher lumi- 
nosity (Compton cooling dominated) and lower luminosity (radia- 
tive cooling dominated) regime is possible. 

We identify these higher and lower luminosity regimes, re- 
spectively, with 'normal' luminosity levels and the occasional 'off' 
states as observed in some X-ray pulsars as e.g. Vela X-l, GX 301- 
2 and 4U1907+05, which are discussed in Section 3. 

In Section 4 we show that the transition from the higher lu- 
minosity to lower luminosity regime may be related to a sudden 
decrease in X-ray photon energy density in the equatorial region 
of the magnetosphere, which is most favoura ble for the plasma t o 
enter due to the Rayleigh-Taylor instability iArons & Lealll97r3) . 
Such a decrease may be the result of a change in the X-ray beam 
pattern from the accretion column (or 'mound' above the polar 
cap) when the X-ray luminosity drops below some critical value 
Lf ~ 3 x 10 35 erg/s determined by the opacity of the column rela- 
tive to Thomson scattering of X-ray photons. Below this luminos- 
ity most of the X-ray emission escapes in a pencil beam, so most 
of the X-ray photons illuminate the magnetospheric cusp region, 
which is stable for plasma entering the magnetosphere. However, 
plasma cooling continues on the longer radiative cooling time scale, 
which is determined only by the density and temperature above 
the magnetosphere, and the source switches into the lower lumi- 
nosity regime. Oppositely, an increase in photon energy density in 
the equatorial magnetospheric region can return the source to the 
higher luminosity regime. 



2 TWO REGIMES OF SUBSONIC ACCRETION FROM A 
QUASI SPHERICAL SHELL 

Consider a hot quasi-spherical shell formed around a slowly rotat- 
ing neutron star magnetosphere, in which accretion proceeds sub- 
sonically (see Paper I for details). The shell can exist as long as 
the X-ray luminosity is less than ~ 4 x 1 36 erg/s, above whi ch 
supersonic (Bondi) accretion is more likely teurnard et aD 19831) . 

To enter the magnetosphere, the plasma in the shell must cool 
down from a high (almost virial) temperature T de termined by hy- 
drosta tic equilibrium [Eq (4) in paper I] to T cr felsner & Lambl 
fl977h 
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Here K is the universal gas constant, fi„, ss 0.6 is the molecular 
weight, G is the Newtonian gravitational constant, M is the neutron 
star mass, k is the local curvature of the magnetosphere, x is the 
angle between the outer normal and the radius-vector at any given 



point at the Alfven surface, and the contribution of turbulent pulsa- 
tions in the plasma to the total pressure is taken into account by the 
factor ( 1 +ym 2 ) (where m, is the turbulent Mach number y = C P /Cy 
is the ratio of specific heat capacities). 

As was shown in Paper I, a transition zone above the Alfven 
surface with radius R A is formed inside which the plasma cools 
down. The effective gravitational acceleration in this zone is 
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and the mean radial velocity of plasma settling is u R = 
/(«) V2GM/ 'R A . The dimensionless settling velocity < f(u) < 1 
is determined by the specific plasma cooling mechanism in this 
zone and is constant through the shell. Together with the density 
of matter near the magnetospheric boundary p(R A ) it determines 
the magnetosphere mass loading rate through the mass continuity 
equation: 



M = 4jxR 2 A p(R A )f(u) yj2GMIR A 



(3) 



This plasma eventually reaches the neutron star surface and pro- 
duces an X-ray luminosity L x = 0. lMc 2 . Below we shall normal- 
ize the mass accretion rate through the magnetosphere as well as 
the X-ray luminosity to the fiducial values M\(, = M/10 16 g/s and 
L36 = Lj/10 36 erg/s, respectively. 

The stationary settling velocity is determined by the plasma 
cooling time t coo \ in the transition zone (see Paper I): 
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where t ff = R 312 / V2GM is the characteristic free-fall time from 
radius R. The angle x is determined by the shape of the magneto- 
sphere, and for the magnetospheric boundary parametrized in the 
form ~ cos A" (where A is the angle counted from the magneto- 
sp heric equator) t an y = n tan A. For example, in model calculations 
by I Arons & Le a ( 1976) n ^ 0.27 in the near-equatorial zone. We 
see that cosx ^ 1 up to A ~ n/2, so below (as in Paper I) we shall 
omit this factor. 



2.1 The Compton cooling regime 

As explained in detail in Paper I (Appendix C and D), in sub- 
sonic quasi-static shells above slowly rotating neutron star magne- 
tospheres such as considered here, the adiabaticity of the accreting 
matter is broken due to turbulent heating and Compton cooling. 
X-ray photons generated near the neutron star surface tend to cool 
down the matter in the shell via Compton scattering as long as the 
plasma temperature T > T x , where T x is the characteristic radia- 
tion temperature determined by the spectral energy distribution of 
the X-ray radiation. For typical X-ray pulsars T x ~ 3 - 5 keV. Cool- 
ing of the plasma at the base of the shell decreases the tempreature 
gradient and hampers convective motions. Addtional heating due to 
sheared convective motions is insignificnt (see Appendix C of Pa- 
per I). Therefore, the temperature in the shell changes with radius 
almost adiabatically KT ~ (2/5)GM/R, and the distance R x within 
which the plasma cools down by Compton scattering is 



R x » 10 u, cm 



V3keV/ 



(5) 



which is much larger than the characteristic Alfven radius R A = 
10 9 cm. 

The Compton cooling time is inversely proportional to the 
photon energy density, 
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t c ~R 2 IL x , 



and near the Alfven surface we find 



t c ~ io[s; 
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This estimate assumes spherical symmetry. Clearly, for the exact 
radiation density the shape of the X-ray emission produced in the 
accretion column near the neutron star surface is important (see 
Sections 3 and 4), but still L x ~ M, Therefore, roughly, f(u)c ~ 
M 1/3 , or, more precisely, taking into account the dependence of R A 
on M in this regime (see Paper I) 
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we obtain: 
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Here yu 30 = yu/10 30 G cm 3 is the neutron star dipole magnetic 
moment. As shown in Paper I, in a spherically symmetric ac- 
cretion flow with turbulence heating and Compton cooling, when 
f(u)c ~ 0.5 (which corresponds to L x ~ 4 x 10 36 erg/s) the sonic 
point emerges above the magnetosphere, accretion becomes super- 
sonic and is most lik ely described by the physical model studied by 
lBumardetalUl983h . 

In the Compton cooling regime a significant X-ray flux vari- 
ability is expected. Suppose that the energy density of X-ray pho- 
tons in the cooling region increases. Then the settling velocity 
f(u)c ~ M >/3 increases as well, leading to further increase in the 
radiation energy density. Clearly, this is a non-stationary situation. 
The maximum accretion rate here will be determined by the ability 
of the magnetosphere to engulf the total amount of mass within the 
cooling region R < R x over the free-fall time tff(R x ): 
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Taking into account the hydrostatic density distribution p(R) 
R~ 3/2 , we find 



Mmax ~ 4nR 2 A p(R A ) ^2GM/R A • | fl - 
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Eliminating density via the mass continuity equation expressed for 
the mean mass accretion rate (M), we obtain 
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(12) 



(M c ) f(u)c 

We stress that this variability occurs on the time scale ~ tff(R x ) 
(about a few hundreds of seconds) and will be present even if the 
density near the magnetosphere is constant. Clearly, variable exter- 
nal conditions would add to the X-ray flux variability on its own 
time-scale. Also note that this instability does not destruct the shell 
as a whole. The reason for this is that above the radius R ~ R x 
Compton heating is effective, so an increase in the X-ray luminos- 
ity would tend to prohibit accretion from the upper layers of the 
shell. 



2.2 The radiative cooling regime 

In the absence of a dense photon field, at the characteristic tem- 
peratures near the magnetosphere T ~ 50-keV and higher, plasma 
cooling is essentially due to radiative losses (bremsstrahlung), and 



the plasma cooling time is r ra j ~ Vt lp- Making use of the con- 
tinuity equation l[3} and the temperature distribution in the shell 
T ~ 1 /R, we obtain 

t mi ~ RM-'fiu). (13) 

Note that, unlike the Compton cooling time J6j, the radiative cool- 
ing time is actually independent of M (remember that M ~ f(u) 
in the subsonic accretion regime!). Numerically, near the magneto- 
sphere we have 
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Following the method described in Section 3 of Paper I, we 
find the mean radial velocity of matter entering the neutron star 
magnetosphere in the near-equatorial region, 
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similar to the expression for f(u) in the Compton cooling region 
Eq. ©. Here £ < 1 is a numerical factor describing the radial ex- 
tension of the transition zone (R A . Using the expression for the 
Alfven radius through f(u), we calculate the dimensionless settling 
velocity: 

m md * o.i^^vr (is) 

and the Alfven radius: 



-6/27 16/27 
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(in the numerical estimates we assume a monoatomic gas with adi- 
abatic index y = 5/3 and turbulent Mach number m, = 1 in the 
shell). The obtained expression for the dimensionless settling ve- 
locity of matter Eq. d 1 6b in the radative cooling regime clearly 
shows that here accretion proceeds much less effectively than in 
the Compton cooling regime (cf. with Eq. 

Unlike in the Compton cooling regime, in the radiative cool- 
ing regime there is no instability leading to an increase of the mass 
accretion rate as the luminosity increases (due to the long charac- 
teristic cooling time), and accretion here is therefore expected to 
occur more quietly under the same external conditions. 

In Fig. Q] we summarize the different regimes of quasi- 
spherical accretion. 



3 'OFF' STATES IN X-RAY PULSARS 

So called 'off' states have been observed in se veral slowly rotat- 
ing low luminosity pulsars such as Vela X-l (Inou eetalJll984l: 
iKrevkenbohm et alJl999Tl2008l : lDoroshenko et alj201ll), GX 301- 
2 dGogiis et alj|201lb and 4U 1907+09 din 't Zand et ail 1 19971: 
ISahiner et alj|2012h ). These states are characterized by a sudden, 
most often without any prior indication, drop in X-ray flux down to 
1 - 10% of normal levels, lasting typically for a few minutes. A few 
examples of lightcurves including off-states are shown in Fig|2] 

It seems to bee fairly well established that the off states can 
not be due simply to increased absorption along the li ne of sight. 



Their occurrence is not c orrelated wi th increased Nh j Fiirst et alj 
l201lllSahineretalJl2012l but see also Kret schmar et al. 1999 for a 



another type of intensity dips in Vela X-l, most probably caused 
by dense blobs in the wind), the tim escale of their onset are too 
short (e.g. IKrevkenbohm et alj [2008) and spectral stu dies show a 
softening of the X-ray spectru m during the off state dGogiis et al.l 
l201ll : [Doroshen ko et alj|201 ll) . contrary to what expected had the 
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REGIMES OF QUASI SPHERICAL ACCRETION 
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Figure 1. Scheme of the different regimes of quasi-spherical accretion including supersonic (Bondi) accretion and subsonic (our model) accretion. 



decreased flux levels been caused by increased absorption. Failure 
by early observations with instruments like RXTEfPCA to detect 
pulsations during the off states have seamed to suggest that the 
sources were instead simply turned off due to the sudden cessa- 
tion of accretion. The popular view is that the cause of this may 
be large density variations in the stellar wind p ossibly combined 
with the onset of the propeller regime (see e.g. lKrevkenbohm et al.l 
l2008t) . 

Recent observations with the more sensitive instruments on- 
board Suzaku of Vela X- 1 ( Doroshenk o et al.l201 ll) . however, show 
that although dropping in luminosity by a factor of ca 20 the source 
is clearly detected with a pulse period equal to that observed at nor- 
mal flux levels. This suggests that rather than cessation of accretion, 
the off-states may be better explained by a transition to a different, 
less effective, accretion regime. We suggest that the onset of the off 
state in these sources marks a transition from the Compton cooling 
dominated to the radiative cooling dominated regime as described 
above. 



4 TRANSITIONS BETWEEN THE COMPTON AND 
RADIATIVE COOLING REGIMES - A CHANGE IN 
BEAM PATTERN? 

In this section we will discuss how transitions between the two 
regimes of plasma entering the magnetosphere may be triggered. 

A decrease in the X-ray photon energy density in the transi- 
tion zone decreases the Compton cooling efficiency, but the Comp- 
ton cooling time remains much shorter than the radiative cooling 
time down to very small luminosities (see Eq. and Eq. d!4t). 
Therefore, in the spherically symmetric case, a transition between 



the two regimes would require an almost complete switch-off of 
the Compton cooling in the equatorial magnetospheric region. In 
the more realistic non-spherical case, the Compton cooling time 
can become comparable to the radiation cooling time when the X- 
ray beam pattern changes with decreasing X-ray luminosity from 
a fan beam to a pencil beam, and the equatorial X-ray flux is re- 
duced by a factor of a few. Additionally, hardening of the pulsed 
X-ray flux with decreasing X-ra y luminosity, which i s observed 
in low-luminosity X-ray pulsars ( Kl ochkov et al.ll20llh . increases 
T x and decreases the specific Compton cooling rate of the plasma 
oc (T — T x )/tc, thus making Compton cooling less efficient. Such 
transitions have been observ ed in transient X-ray pulsars (see, e.g., 
jParmar. White & Stella 1989)). The radiation density in the pencil 
beam cools down the plasma predominantly in the magnetospheric 
cusp region, but be cause of the stronger magnetic line curvature 
( Arons & Lea 1976) the plasma entry rate through the cusp will 
be insignificant. Still, the plasma continues to enter the magneto- 
sphere via instabilities in the equatorial zone, but at a lower rate 
determined by the longer radiative cooling timescale. 

The mass accretion rate in the radiative cooling regime will 
be determined by the plasma density by the time Compton cooling 
switches off in the magnetospheric equator region. This occurs at 
some X-ray luminosity L x < L„, 

In the case of a strong neutron star magnetic field (> 
10 12 G) most of the thermal X-ray photons produced in the 
energy release zone are produced with the ordinary (O) polar- 
ization mode, and the number of extraordinary (X) photons is 
small. Depending on the plasma density, the vacuum polariza- 
tion effects leading to conversion of O-photons into X-photons 
can be important for photons with energies between fiu> v » 
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Figure 2. Example lightcurves of Vela X-l, 4U 1907+09 and GX 301-2 (HEASARCH archive data). Shaded areas mark observed off-states. Note that the 
timescale is different in the different panels. 



12(« e /10 22 cm" 3 ) 1/2 (B/10 12 Gr 1 keV and the cyclotron resonance 



energy E c = fiu) « 1 1 6(5/10 12 G) keV j Ventura et al J 1 19791 ; 
iMeszaros & Nagell 1 19851 ; lHarding & Lap l2006h . From equating 
tuo v = tuo c , we find the critical electron number density 



7.6 x io- [cm-] yy 



(18) 



below which vacuum polarization effects on the photon mode prop- 
agation are significant. Since n e ~ M/A ~ MR A (here A ~ 
nRj /s 9 2 c ~ R 3 NS /R A is the effective area of the NS surface onto which 
accretion proceeds) increases with mass accretion rate, in bright X- 
ray pulsars with L x > 10 37 erg/s vacuum polarization effects are not 
expected to play a significant role, and the X-ray emission beam 
should consist mostly of O-photons. 

O-photons with energies below the cyclotron resonance en- 
ergy E c propagate along the magnetic field with a scattering cross- 
section cry « cr T sin 2 1? where & is th e angle between the photon 
wave vector and the magnetic field ( Harding & Lai 2006). They 



form a pencil beam with a characteristic opening angle 9 P ~ 
1/ Vtt, where t t is the optical dep th in the energy release zone 
(Baskd ll976l ; [Dplginov et al,lll979l) . These photons exert a low 
force on the accreting matter, and if the number of hard photons 
with energies above E c (where cry » cr T ) is small, no high accre- 
tion column will be formed. Apparently, this is the case in some 
luminous X-ray pulsars (Cen X-3, Her X-l) with L x ~ 10 37 erg/s. 

At lower densities (corresponding to the low mass accretion 
rates and low X-ray luminosities we are considering here) or higher 
magnetic fields, vacuum polarization is significant, and O-photons 
can be converted into X-photons. The scattering cross-section of 
such photons is cr ± = a- T (E/E c ) 2 and is independent of the angle 
&. Therefore, they are expected to produce a more spherically sym- 
metric (but still not fan-like) beam. Again, no high accretion col- 
umn is expected to form if the number of hard photons with E > E c 
is small. 

Additionally, photons must be scattered in the accretion flow 
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above the polar cap region. This scattering does not significantly 
affect the pencil beam formed by O-photons in bright pulsars be- 
cause of the sin 2 & dependence of the scattering cross-section. In 
low-luminosity pulsars, in which vacuum polarization photon mode 
conversion occurs, Thomson scattering would tend to form a fan 
beam at all photon energies. Therefore, the transition from fan to 
pencil beam in low-luminosity X-ray pulsars without high columns 
does not occur until the optical depth in the accretion flow above 
the polar cap becomes less than one. The optical depth in the ac- 
cretion flow in the direction normal to the neutron star s urface from 
the ra dial distance r$ = r/10 6 cm is estimated to be dLamb et"al] 
Il973l) 

(here the neutron star mass is assumed to be 1.5 M Q and the NS 
radius R NS = 10 6 cm). Taking into account the dependence of the 
Alfven radius on M and n ([8), we see from this estimate that the X- 
ray diagram change is expected to occur at r,. < 1, corresponding 
to an X-ray luminosity of 

L t ~ 3 x 10 35 [erg/s]^ 3/I °. (20) 

The decrease in radiation energy density in the magneto- 
spheric equatorial zone due to the X-ray beam pattern change leads 
to an increase of the Compton cooling time and hence triggers a 
transition to the lower luminosity regime and the source enters the 
'off' state. From the mass continuity equation we then find the lu- 
minosity ratio: 

L x ,rad _ f(u)md I k \ 3 ^1) 
/(")C Vrarf/ 

Substituting expressions for /(«)'s Eq. (9J and Eq. l |16l > taken at 
L x = Lf, we find the X-ray luminosity in the 'off' state: 

L„ad« 10 35 [erg/s]^ 33 . (22) 

We stress here that this X-ray luminosity is derived for the 
case of a shell around the neutron star magnetosphere in the ra- 
diation cooling regime. Lower X-ray luminosities can be realized 
only if the density of matter near the magnetosphere p(Ra) (which 
is determined by the density of gravitationally captured stellar wind 
behind the shock at the Bondi radius p(Rg)) turns out to be smaller 
than the value that provides the minimum X-ray luminosity ~ L\ 
of the source, which is required for effective Compton cooling to 
operate in the equatorial region of the NS magnetosphere. 

The return from radiative cooling dominated accretion back to 
the Compton cooling dominated regime can take place, for exam- 
ple, due to a density increase above the magnetosphere, leading to 
an increase of the mass accretion rate. In turn, this leads to growth 
of the vertical optical depth of the accretion column, disappear- 
ance of the beam and enhancement of the lateral X-ray emission. 
The radiative energy density in the equatorial magnetospheric re- 
gion strongly increases, Compton cooling resumes, and the source 
comes back to normal luminosity levels. 

4.1 The changing beam pattern in Vela X-l 

The idea that the transition between the two regimes may be trig- 
gered by a change in the X-ray beam pattern is supported by ob- 
servations of Vela X-l. In Fig [3] we plot the Suzaku pulse pro- 
fil es of Vela X-l in differe nt energy bands from the observation 
bv lDoroshenko et ail d201ll) . The top 4 panels show the pulse pro- 
files at normal luminosity levels and the bottom 4 panels show the 



pulse profiles at different energies within an off-state. The observed 
change in phase of the 20-60 keV profile in the off-state (at X-ray 
luminosity ~ 2.4 x 10 35 erg/s) suggest a disappearance of the fan 
beam at hard X-ray energies upon the source entering this state. 
Instead, a pencil beam formed by thermal O-photons from the two 
magnetic poles of the neutron star is clearly visible. This pencil 
beam is schematically shown as a dashed region superimposed on 
the other profiles in all panels . 

It is of importance here that a cyclotron resonance feature in 
Vela X-l is observed around 20 keV, and that above this energy 
the scattering cross-section of both O- and X-photons is essentially 
equal to the Thomson value. This is why the 20-60 keV pencil beam 
profile seen in the 'off' state (where the X-ray luminosity L < L t ) is 
strongly scattered by the Thomson-thick accretion flow above the 
polar caps to form a fan beam at normal luminosity levels. Around 
20 keV resonance scattering in the cyclotron line above the polar 
cap is significant, as suggested by the 12-20 keV profile in the 'off' 
state (panel number 6 from the top in Fig.[3]l. At lower energies the 
absorption in the pencil beam dominates in the 'off' state as well as 
at normal luminosity levels. 



5 DISCUSSION 

We have considered two possible regimes of plasma entering the 
magnetosphere - Compton-cooling dominated and radiative cool- 
ing dominated subsonic accretion, corresponding to X-ray lumi- 
nosities differing by more than an order of magnitude. It is essential 
to realize that the two regimes can be realized for the same density 
of plasma around the magnetosphere. The density p(Ra) at the bot- 
tom of the quasi-static shell near the magnetosphere depends on 
the density behind the external shock in the wind at the Bondi cap- 
ture radius R B as p(Ra) ~ p(Rb)(Rb/Ra) 3/2 ~ P« v ^ 3 > i- e - it is vel 7 
sensitive to the stellar wind density and, especially, to its velocity. 

On the contrary, the accretion rate onto the neutron star is in 
our model determined by the ability of the plasma to enter the 
magnetosphere, i.e. by the plasma settling velocity. If the radia- 
tion energy density is high enough and the Compton cooling time 
is shorter than the radiative cooling time of the plasma, the mean 
settling plasma velocity f(u) is about 0.3 times the free-fall veloc- 
ity. If the radiation density decreases and the radiative cooling time 
becomes shorter than the Compton cooling time, the mean settling 
velocity of the plasma decreases by a factor of 3, and the mass en- 
try rate onto the magnetosphere drops by a factor of ~ 4 relative to 
that due to Compton cooling at a given density p(Ra)- 

A likely mechanism responsible for the triggering of a tran- 
sition between the two regimes is the change of the X-ray emis- 
sion diagram from the accretion column from being essentially a 
fan beam at high luminosities (large optical thickness of the ac- 
cretion column) to a pencil beam (small vertical optical thickness 
of the accretion column) when the mass accretion rate decreases 
below ~ 3 x 10 15 g/s, corresponding to an X-ray luminosity of 
L| = 3 x 10 35 erg/s. For the near-magnetosphere density corre- 
sponding to this X-ray luminosity the stationary mass entry rate 
due to radiative cooling is about 10 15 g/s. Lower mass accretion 
rates are possible for lower external wind densities. 

Clearly, the complicated picture of accretion onto magnetized 
rotating neutron stars is far from being complete. In our model 
we have taken into account turbulent heating as well as Comp- 
ton and radiative cooling of the plasma near the magnetosphere, 
but ignored the possible effects of magnetic fields frozen into the 
plasma, which should be studied separately (see, e.g., resent stud- 
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Figure 3. Pulse profiles of Vela X-l as observed by Suzaku (in normalized counts from Doroshenko e t alj20l"lh at normal luminosity levels (4 upper panels) 
and in an 'off' state (4 lower panels). The dashed region schematically shows the pencil X-ray beam seen from the magnetic poles of the neutron star in the 
'off' state at 20-60 keV (panel 5 from the top) superimposed on the other profiles. 



ies bv llkhsanov et al.l ([2012)). Magnetic field reconnection, on the 
one hand, serves as a source of additional heating of the plasma, 
but on the other hand, may f acilitate plasma entering the magneto- 
sphere ( Eisn er & Lambl 19841) . The magnetic reconnection near the 
magnetosphere may be responsible for occasional transitions into 
the 'strong coupling regime', in which drastic jumps of the neutron 
star spin period may occur without any change in X-ray luminosity, 
as discussed in Paper I. This issue should be further investigated. 



6 CONCLUSIONS 

To conclude, quasi-spherical accretion from a stellar wind onto a 
slowly rotating magnetized neutron star may proceed in different 
ways. In high-luminosity sources with L x > 4 x 10 36 erg/s the mat- 
ter behind the bow shock at the Bondi radius cools down rapidly 
and falls freely towards the magnetosphere, forming a shock above 
the Alfven surface as described by models of Bondi- or supersonic 
accretion. 

At lower X-ray luminosities a quasi-static atmosphere is 
bound to be formed above the neutron star magnetosphere, and ac- 
cretion proceeds subsonically with an accretion rate determined by 
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the ability of the plasma to enter the magnetosphere via instabili- 
ties. A model for such subsonic quasi spherical accretion was pre- 
sented in Paper I. 

At luminosities above L t ~ 5 x 10 35 erg/s, the plasma cools 
down via Compton processes and enter the magnetosphere in the 
equatorial regions, most favorable for the Rayleigh-Taylor insta- 
bility to develop. At X-ray luminosities below L t ~ 5 x 10 35 erg/s, 
however, the radiation energy density in the magnetosphere equator 
may be significantly reduced due to a change in the X-ray emission 
pattern from the accretion column. At L > L t the optical depth of 
the accretion column is larger than one and radiation is scattered 
mostly in the lateral direction forming a fan beam, at L x < L t the 
vertical Thomson optical depth in the accretion column becomes 
smaller than one, and a pencil-beam is formed. The pencil beam 
illuminates the magnetosphere cusp region, where plasma entry is 
hampered by large curvature of the magnetic field lines. Still, the 
plasma continues to find its way onto the magnetic field lines in 
the equatorial region due to radiative plasma cooling, which is de- 
termined by the plasma density independently of the radiation en- 
ergy density. The mass accretion rate due to radiative cooling is by 
several times smaller than due to Compton cooling with the same 
plasma density above the Alfven surface. 

We identify the two subsonic regimes; the Compton cooling 
dominated and the radiative cooling dominated, respectively, with 
observed 'normal' luminosity levels and the so called 'off' states 
in some slowly rotating X-ray pulsars, e.g. Vela X-l, GX 301-2, 
4U1907+05. Our proposed scenario is supported by the observed 
change in the hard X-ray pulse profile in the off state of Vela X- 1 
and could be further checked against observations of the behaviour 
of hard X-ray pulse profiles during a transition into or out of an 
off-state. 

It also can not be excluded that the phenomenon of Super- 
giant Fast X-ray Transients (SFXTs) (see ISidolil I20TTI for a re- 
cent summary and review) can similarly be related to transitions 
between different regimes of plasma cooling in a quasi-spherical 
shell around a slowly rotating magnetized neutron star. The qui- 
escent states of SFXTs with stable low-luminosity accretion with 
Lrad ~ 10 34 erg/s may be controlled by thermal plasma cooling, 
while the unstable X-ray flares may be triggered by a density in- 
crease above the magnetosphere leading to an increase in the lat- 
eral X-ray emission from the accretion column and a transition to 
the Compton-cooling dominated regime. Our model predicts the 
corresponding change in hard X-ray pulse profile to occur during 
the transition from the quiescent state to the flaring state, which can 
be checked by dedicated observations. 
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